Predictions for Neutron Scattering and Photoemission Experiments on CuGeOs 
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Applying numerical techniques to a model recently proposed for the one dimensional (ID) spin- 
Peierls compound CuGeC>3, we calculate dynamical properties that can be directly compared with 
inelastic neutron scattering (INS) data and angle-resolved photoemission (ARPES) experiments. 
The momentum and energy dependence of the dynamical structure factor S(q, to) are discussed, as 
well as the static structure factor S(q). The spectral function A(q, to) calculated from the one particle 
Green's function at half-filling is shown at several values of the hole hopping amplitude t. The results 
have some unique features characteristic of one dimensional systems including small weight near the 
Fermi energy. The presence of "shadow bands" induced by short distance antiferromagnetic (AF) 
correlations is predicted to appear in ARPES experiments for CuGeC>3 and also for Sr2CuC>3. 
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The behavior of low dimensional magnetic materials 

^Continues to attract much experimental and theoretical 

^interest. A variety of exotic ground states can be ob- 
"^nserved in different compounds depending on anisotropies, 
■^^the spin of the magnetic ions, and the coupling to the lat- 
i-Htice degrees of freedom. Recently, CuCeC^ was found to 
be the first inorganic system that shows a spin-Peierls 
OO(SP) phase. [1] With decreasing temperature a regime 
"^Fdescribed by uniform S=l/2 Heisenberg chains under- 
^^goes a phase transition at Tsp = 14K into a dimerized 
(< — ^system due to the coupling of the spin- 1/2 Cu 2+ -ions to 
ly-jthe three-dimensional lattice phonons. At Tsp a struc- 
^\tural transition takes place where alternate atoms are 
"^displaced in opposite directions, and an energy gap for 
agnetic spin triplet excitations appears. The results of 



QHase, Terasaki and Uchinokura [1] have generated much 
..^excitement, and the presence of a spin-gap in CuCeC^ 
(— ih as been confirmed by several groups. [2,3] 
Q CuGeOs has a c-direction lattice constant of 2.941A, 
Omuch smaller than the a and b lattice parameters. Al- 
though the c-axis lattice constant corresponds to the 
Cu 2+ — Cu 2+ distance, the actual links between cop- 
pers are provided by the edge sharing of CuOe octa- 
hedra. The angle defined by each Cu-O-Cu bond is 
~ 86° (Ref. [4]). This induces an exchange coupling Ji 
in the effective nearest-neighbors (n.n.) S=l/2 Heisen- 
berg model for the Cu-0 chains which is much smaller 
than those found in the two dimensional (2D) cuprate 
superconductors. The measured susceptibility, above 
14K while showing characteristics of ID antiferromag- 
nets, is not quantitatively reproduced by calculations cor- 
responding to a n.n. S=l/2 Heisenberg model. Actually, 
Lorenzo et al. [5] based on a neutron scattering study 
suggested that a next-nearest-neighbor (n.n.n.) coupling 
J2 may arise from the Cu-O-O-Cu path. [5] Riera and 
Dobry [6] arrived to the same conclusion studying nu- 
merically the temperature dependence of x- They re- 
ported a substantial ratio J2/J1 = 0.36. In this regime 



both phonons and frustration contribute to the open- 
ing of the spin gap. Castilla et al. [7] also remarked 
the relevance of n.n.n. terms reporting a smaller ratio 
J2/J1 = 0.24, with Ji = 150K, which is right below the 
ratio where the effects of frustration open a spin gap in 
the model in the absence of phonons. Actually, exact 
diagonalization (ED) techniques locate the critical ratio 
at J2/J1 |c ~ 0.2412 ± 0.0001. [8,7] While the microscopic 
origin of J2 is not quite clear, we include it in our studies 
simply on phenomenological grounds to model accurately 
CuGeOs. Since some experimental evidence shows that 
the SP transition is exclusively driven by phonons, [3] in 
our study we use Castilla et al.'s estimation of J2. How- 
ever, note that working so close to J2/J1IC likely induces 
strong dimerizing fluctuations of non-phononic origin in 
the ground state. 

Based on this discussion, the model studied here cor- 
responds to a Heisenberg model with n.n. and n.n.n. in- 
teractions and a term that dimerizes the lattice to mimic 
the effects of phonons. The Hamiltonian is 

H = Ji^Si-Sj+J 2 J2 S 1 -S m + ^(-l) 1 S 1 -S J , 

(y> «im» (ij) 

(1) 

where i,j,m denote sites of a ID chain with N sites, S; 
are S=l/2 spins, and (ij) (((im))) corresponds to n.n. 
(n.n.n.) spins along the chain. The static properties of 
the special case J 2 = have been widely studied in the 
literature, [9] while the frustrated Ji — J2 model has also 
received much attention. [8] However, dynamical studies 
with controlled approximations have not been previously 
reported to the best of our knowledge, and here we pro- 
vide this information using ED techniques. [10] 

Inelastic neutron scattering experiments have been re- 
cently performed on CuGeOs single crystals. [14,15] The 
opening of a gap at q = was observed in the dispersion 
of the magnetic excitations. An ED study by Castilla 
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et al. [7] showed that the INS triplet dispersion is re- 
produced using Eq.(l) if J 2 /Ji = 0.24 and 6 = 0.03 Ji. 
We adopt these parameters in the present paper. Thus 
far the dynamical information contained in the energy 
dependence of S(q, w) has not been analyzed in the INS 
literature for CuGeC>3. The reason is that the experimen- 
tal data is usually contaminated by spurious effects, and 
a complicated convolution between the actual data and 
the instrumental resolution is needed to extract S(q, w). 
To guide this experimental effort it is useful to have a 
theoretical prediction for S(q, w). The purpose of this 
paper is to provide for the first time spin dynamical in- 
formation obtained from model Eq.(l) to compare the- 
ory with experiments in CuGeC>3. Here we also calcu- 
late the spectral function A(q, w) which can be directly 
compared with ARPES experiments. In the context of 
high temperature superconductors, recent ARPES data 
for Sr2Cu02Cl2 at half-filling has provided important 
constraints on theoretical models for the cuprates. [11] 
A similar study in CuGeOs would likely prove as useful. 
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FIG. 1. (a) S zz (q,i 
with a broadening e = 



40.0 



0.0 



(b) 



0.0 0.5 1.0 

q/n 



continuum 



0.0 0.2 0.4 0.6 0.8 1.0 

•) obtained on a chain of N = 20 sites 
0.3Ji at several momenta; (b) Position 
of the first pole in Fig. la for each q, defining the spin triplet 
dispersion (diamonds). The circles denote results from Refs. 
[14,15]. The dashed line near these points guides the eye , 
while the dashed line at higher energies signals the position 
of the pole with the highest energy for each q. In the inset, 
the intensity h of the highest peak in Fig. la, normalized with 
respect to q = 7r, is shown as a function of momentum, using 
J2/J1 = 0.24, S/h = 0.03 and Ji = 150K. 

In Fig. la S zz (q,w) is shown at several momenta. The 
position of the dominant peak defines the energy of the 
first triplet excitation for each q. [12] Note that this 
peak is not isolated but forms part of a broad contin- 
uum. The intensity is the highest at q = 7r as expected 
due to the power-law decaying AF spin-spin correlations 
in the ground state. The intensity of the dominant peak 
for each momentum, normalized with respect to q = 7r, is 
shown in the inset of Fig. lb. This ratio decays away from 
q = 7r more slowly than in the case of a S=l Heisenberg 



chain [13] which has a spin gap caused by a disordered 
but translationally invariant ground state. In Fig. lb, the 
first peak in the spectrum is shown as a function of mo- 
mentum and compared with experiments. [14,15]. As 
previously remarked by Castilla et al. [7] the agreement 
between theory and experiment is excellent. From our 
numerical study we can also analyze the position of the 
peak with the highest energy to put bounds on the mul- 
tiparticle continuum. [16] The result is shown in Fig. lb, 
which in principle could be compared directly with INS 
results. However, note that the rapidly decaying inten- 
sity of the spectral weight makes such comparison diffi- 
cult. If the resolution of the INS experiment only allows 
for the observation of intensities within, e.g., 10% or 20% 
of the maximum at q = 7r, then the apparent width of the 
continuum would appear considerably reduced. 
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FIG. 2. S(q), normalized with respect to its value at q = 7r, 
as a function of q. The convention for the symbols is shown 
in the inset. We present results for the model Eq.(l), for the 
(gapless) n.n. spin-1/2 Heisenberg chain, and for the spin-1 
Heisenberg chain which has a disordered ground state and a 
finite spin gap. 

In Fig. 2, the static structure factor S(q) ( = f du> 
S(q, w)) is presented. As shown in the figure, finite size 
effects are small and we consider our results representa- 
tive of the bulk limit. As expected, the large AF corre- 
lations in the model cause S(q) to peak at q = 7r, slowly 
decreasing away from it as the momentum is varied. The 
normalized results for Eq.(l), with the parameters intro- 
duced in Ref. [7], can be clearly distinguished from those 
obtained at J2 = 0, and thus INS measurements of S(q) 
can be used to confirm the presence of important n.n.n. 
couplings in CuGeOs (8 is so small that does not influ- 
ence on S(q) as much as J2). For completeness, results 
for the S=l chain are also shown. 

We have also calculated the spectral function A(q, w) 
which corresponds to the creation of holes by the removal 
of spins on the chains in the sudden approximation. As 
Hamiltonian for this calculation we use the t — J model 
generalization of Eq.(l) i.e. we add a hopping term 
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Ht = -t$>kc ji<T +h.c.), (2) 

where c are hole operators and the rest of the notation 
is standard. [10] In the absence of experiments testing 
the properties of carriers in CuGeOs, it is difficult to 
predict the value of the hopping amplitude t. Thus, in 
the results quoted below three different amplitudes are 
used such that more information is available to compare 
with experiments once ARPES data becomes available. 
However, if the couplings Ji and J 2 are caused by an ex- 
change mechanism it is likely that the amplitude t would 
be larger than both. For simplicity, we have neglected 
the influence of the dimerization on t which in princi- 
ple should also be modulated as the exchange Ji. Since 
the dispersive features observed in our calculations have 
bandwidths of order Ji (>> 6 ), the influence of 6 on the 
ARPES data is expected to be small. 
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FIG. 3. Spectral function A(q, to) on a chain of 20 sites, 
at several values of t/Ji. The broadening is e = 0.25Ji, and 
the Fermi energy is arbitrarily located at to = 0. The other 
couplings are as in Fig.l. 

In Fig. 3, A(q, w) is shown for different momenta and 
three couplings t/Ji, with Ji, J2 and 6 defined as before. 
Note the presence of a dominant peak at low binding en- 
ergy which moves towards the Fermi energy (at u> = in 
the figure) as the momentum is changed from q = to 
q = 7r/2 for all couplings. The lowest binding energy is 
precisely reached at q = 7r/2. Based on previous experi- 
ence with the 2D copper oxide compounds, [10] this large 
peak may correspond to a hole immersed in a locally or- 
dered AF spin background which increases its effective 
mass. An interesting feature of our results is the pres- 
ence of substantial weight at large u> for all q and cou- 
plings. At q = 7r/2, the shape of A(q, w) at the top and 
bottom of the band seems similar i.e. a symmetry line 
seems to run through the middle of the total bandwidth. 
This feature also appears in the 2D t — J model but only 
at very small J/t. [17] This is another effect predicted 
by theoretical calculations that an ARPES experiment 



could analyze, although we are aware that in such ex- 
periments it is difficult to investigate structure at large 
binding energies due to the presence of core backgrounds 
in the data that grow with binding energy, as well as the 
influence of other bands not taken into account in simple 
one band models like the one used in the present paper. 

A curious effect obvious to the eye in Fig. 3 is the pres- 
ence of a robust peak in A(q, u>) for momenta larger than 
q = 7r/2. This behavior is quite different from what is 
observed in normal metals, and here we argue that it 
is caused by the strong AF correlations in the ground 
state. Recently, this effect, which is usually referred to 
as "shadow bands" , has received much attention in the 
context of the 2D cuprates after ARPES experiments ad- 
dressed their existence in Bi2212 at optimal doping. [18] 
The origin of the shadow bands is simple: in the pres- 
ence of AF order the unit cell doubles its size due to the 
nonzero staggered order parameter, reducing the Bril- 
louin zone in half. This induces an extra symmetry in 
the spectrum since now, e.g., on a square lattice mo- 
menta along the diagonal from (0,0) to (7r/2,7r/2) be- 
come equivalent to those from (7r/2,7r/2) to (7r,7r). [19] 
ARPES techniques [11] have shown that antiferromag- 
netically induced bands can be clearly seen in the insula- 
tor Sr2Cu02Cl2- If AF order is replaced by a finite but 
robust AF correlation, the effects of shadow bands are 
expected to be still visible in ARPES data. The present 
calculations have shown that shadow bands can be ro- 
bust even in ID systems, where we are used to the con- 
cept that there are no ordered phases at zero tempera- 
ture. Our prediction can be extended to materials like 
Sr2CuC>3 where the n.n. Heisenberg model accurately de- 
scribes its properties. [20] The existence of shadow bands 
in ID systems should not be surprising since the hole dis- 
persion tested in ARPES experiments is dominated by 
short distance correlations. 

In Fig.4a,b,c, the dispersion of the large features 
(peaks) observed in Fig. 3 are plotted, both at low and 
high binding energy. The approximate symmetry of the 
spectrum with respect to the middle of the total band- 
width is clear. In Fig.4d the weight Z of the lowest en- 
ergy pole appearing in A(7r/2,w) is shown as a function 
of N -1 / 2 as presented in previous studies for ID systems. 
[21] Although Z is very small, it nevertheless seems to 
converge to a finite constant in the bulk even at very 
large ratios t/Ji contrary to what occurs for the pure 
t — J model in ID. [21] Emery [22] recently argued that 
in the present calculation at half — filling with a spin 
gap in the system, the charge degrees of freedom tested in 
A(q, w) indeed should have a one particle Green's func- 
tion with a quasiparticle peak carrying a finite Z (see 
also Ref. [21]). The strong correlations reduce Z to only 
a fraction of its maximum value 1 (as shown in Fig.4d), 
but ARPES experiments should still observe a peak with 
a finite weight in their data. However, as the density 
of carriers becomes finite, the holes will tend to form a 
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Luttinger liquid with Z=0. The crossover between these 
two regimes is interesting both theoretically and experi- 
mentally if CuGeC>3 is doped away from half-filling, and 
it certainly deserves further study. 
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FIG. 4. Dispersion of the poles with the lowest and high- 
est binding energy found in A(q, to) using a 20 site chain at 
different q and for coupling t/Ji equal to (a) 0.5, (b) 1.0 and 
(c) 2.0. The rest of the couplings are as in Fig.l; (d) Spec- 
tral weight Z(tt/2) of the pole with the lowest binding energy 
appearing in our calculation at q = 7r/2 vs N -1 ^ 2 , where N 
is the number of sites. The triangles, diamonds, squares and 
circles correspond to t/Ji equal to 2.0, 8.0, 32.0 and 128.0, 
respectively. Z(?r/2) is normalized such that it is bounded 
between and 1. 

Summarizing, here we presented dynamical properties 
of a model for CuGeC>3. Theoretical predictions are made 
that can be directly compared with experiments. In the 
context of INS, we have calculated the intensity and mo- 
mentum dependence of S(q, u>), as well as the integrated 
static structure factor S(q). For ARPES experiments at 
half-filling, the dispersion of quasiparticle-like features 
was provided. The shape of the spectrum is somewhat 
unusual with substantial weight concentrated at high en- 
ergies. We also predict the presence of "shadow bands" 
in ARPES data induced by AF correlations in the ground 
state of the ID compounds CuGeOs and S^CuOs. 
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